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On the exact gravitational lens equation in spherically symmetric and static 

spacetimes 

Volker PerlidQ 
Institut fur Theoretische Physik. Universitat zu Koln 
Ziilpicher Str. 77, 50937 Koln, Germany Q 

Lensing in a spherically symmetric and static spacetime is considered, based on the lightlike 
geodesic equation without approximations. After fixing two radius values ro and rs, lensing for 
an observation event somewhere at ro and static light sources distributed at rs is coded in a lens 
equation that is explicitly given in terms of integrals over the metric coefficients. The lens equation 
relates two angle variables and can be easily plotted if the metric coefficients have been specified; 
this allows to visualize in a convenient way all relevant lensing properties, giving image positions, 
apparent brightnesses, image distortions, etc. Two examples are treated: Lensing by a Barriola- 
Vilenkin monopole and lensing by an Ellis wormhole. 

PACS numbers: 04.20.-q, 98.62.Sb 



I. INTRODUCTION 

Theoretical work on gravitational lensing is tradition- 
ally done in a quasi-Newtonian approximation formal- 
ism, see, e.g., Schneider, Ehlers and Falco [l| or Petters, 
Levine and Wambsganss which is based, among other 
things, on the approximative assumptions that the grav- 
itational field is weak and that the bending angles are 
small. Under these assumptions, lensing is described in 
terms of a "lens equation" that determines a "lens map" 
from a "deflector plane" to a "source plane" , thereby re- 
lating image positions on the observer's sky to source 
positions. Although for all practical purposes up to now 
this formalism has proven to be very successful, there 
are two motivations for doing gravitational lens theory 
beyond the quasi-Newtonian approximation. First, from 
a methodological point of view it is desirable to inves- 
tigate qualitative features of lensing, such as criteria for 
multiple imaging or for the formation of Einstein rings, in 
a formalism without approximations, as far as possible, 
to be sure that these features are not just reflections of 
the approximations. Second, lensing phenomena where 
strong gravitational fields and large bending angles are 
involved are no longer as far away from observability as 
they have been a few years ago. In particular, the dis- 
covery that there is a black hole at the center of our 
galaxy Q, and probably at the center of most galaxies, 
has brought the matter of lensing in strong gravitational 
fields with large bending angles closer to practical astro- 
physical interest. If a light ray comes sufficiently close to 
a black hole, the bending angle is not small; in principle, 
it may even become arbitrarily large, corresponding to 
the light ray making arbitrarily many turns around the 
black hole. Unboundedly large bending angles also oc- 
cur e.g. with wormholes; the latter are more exotic than 
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black holes, in the sense that up to now there is no clear 
evidence for their existence, but nonetheless considered 
as hypothetical candidates for lensing by many authors. 

If one wants to drop the assumptions of weak fields 
and small angles, gravitational lensing has to be based 
on the lightlike geodesic equation in a general-relativistic 
spacetime, without approximations. In this paper we 
will discuss this issue for the special case of a spheri- 
cally symmetric and static spacetime. In view of applica- 
tions, this includes spherical non-rotating stars and black 
holes, and also more exotic objects such as wormholes 
and monopoles with the desired symmetries. The main 
goal of this paper is to demonstrate that in this case lens- 
ing without approximations can be studied, quite conve- 
niently, in terms of a lens equation that is not less explicit 
than the lens equation of the quasi- Newtonian formalism. 

Lensing without weak-field or small-angle approxima- 
tions was pioneered by Darwin 0,O| and by Atkinson |(|. 
Whereas Darwin's work is restricted to the Schwarzschild 
spacetime throughout, Atkinson derives all relevant for- 
mulas for an unspecified spherically symmetric and static 
spacetime before specializing to the Schwarzschild space- 
time in Schwarzschild and in isotropic coordinates. All 
important features of Schwarzschild lensing are clearly 
explained in both papers. In particular, they discuss the 
occurrence of infinitely many images, corresponding to 
light rays making arbitrarily many turns around the cen- 
ter and coming closer and closer to the light sphere at 
r = 3m. However, they do not derive anything like a 
lens equation. 

The notion of a lens equation without weak-field or 
small-angle approximations was brought forward much 
later by Frittelli and Newman 0. It is based on the 
idea of parametrizing the light cone of an arbitrary ob- 
servation event in a particular way. For a general discus- 
sion of this idea and of the resulting "exact gravitational 
lens map" in arbitrary spacetimes the reader may con- 
sult Ehlers, Frittelli and Newman Q or Per lick Here 
we are interested only in the special case of a spherically 
symmetric and static spacetime. Then the geodesic equa- 



tion is completely integrable and the exact lens equation 
of Frittelli and Newman can be written quite explicitly. 
One can evaluate this equation from the spacetime per- 
spective, as has been demonstrated by Frittelli, Kling 
and Newman ^(J for the case of the Schwarzschild space- 
time, thereby getting a good idea of the geometry of the 
light cone. Here we will use an alternative representa- 
tion, using the symmetry for reducing the dimension of 
the problem. After fixing two radius values ro and rs, 
lcnsing for an observation event somewhere at tq and 
static light sources distributed at r$ is coded in a lens 
equation, explicitly given in terms of integrals over the 
metric coefficients, that relates two angles to each other. 
This representation results in a particularly convenient 
method of visualizing all relevant lensing properties, as 
will be demonstrated with two examples. 

The lens equation discussed in this paper should be 
compared with the lens equation for spherically symmet- 
ric and static spacetimes that was introduced by Virb- 
hadra, Narasimha and Chitre fill and then, in a modified 
form, by Virbhadra and Ellis [12] • The Virbhadra-Ellis 
lens equation has found considerable interest. It was ap- 
plied to the Schwarzschild spacetime [12( and later also to 
other spherically symmetric and static spacetimes, e.g. to 
a boson star by Dabrowski and Schunck Il3l| , to a fermion 
star by Bilic, Nikolic and Viollier to spacetimes 

with naked singularities by Virbhadra and Ellis |l5|. to 
the Reissner- Nordstrom spacetime by Eiroa, Romero and 
Torres and to a Gibbons-Maeda-Garfinkle-Horowitz- 
Strominger black hole by Bhadra 0]. In the last two 
papers, the authors concentrate on light rays that make 
several turns around the center and they use analytical 
methods developed by Bozza The Virbhadra-Ellis 

lens equation takes an intermediary position between the 
exact lens equation and the quasi-Newtonian approxima- 
tion. It makes no assumptions as to the smallness of 
bending angles, but it does make approximative assump- 
tions as to the position of light sources and observer. For 
the Virbhadra-Ellis lens equation to be valid the space- 
time must be asymptotically flat for r — > oo and both 
observer and light sources must be at positions where r is 
large; moreover, one has to restrict to light sources close 
to the radial line opposite to the observer position, i.e., 
to the case that there is only a small misalignment. (The 
question of how one can free oneself from the latter as- 
sumption was addressed by Dabrowski and Schunck [l3| 
and by Bozza [19|.) The lens equation to be discussed in 
the present paper is not restricted to the asymptotically 
flat case, and it makes no restriction on the position of 
light sources or observer. 



II. DERIVATION OF THE LENS EQUATION 

We consider an arbitrary spherically symmetric and 
static spacetime. For our purpose it will be advantageous 
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FIG. 1: The figure shows the plane i? = tt/2. The observer, 
indicated by a dot, is situated at r = ro, f = 0, the light 
sources are distributed at r = rs- The lens equation relates 
the angle which gives the image position on the observer's 
sky to the angle <£> which gives the source position in the 
spacetime. The angle ip indicated in this picture is irrelevant 
for the lens equation but will be needed in Subsection II V Cl 

to write the metric in the form 

g = A{rf ( S{rf dr 2 + R{rf (dtf 2 +sin 2 z9 dip 2 ) - dt 2 ^. 

Here ip and # are the standard coordinates on the sphere, 
t ranges over R and r ranges over an open interval 
] fmin , r max [ where -oo < r min < r max < oo. We as- 
sume that the functions A, S, and R are strictly posi- 
tive and (at least piecewise) differentiable on the inter- 
val ] r min , r max [ . As the lightlike geodesies are not 
affected by the conformal factor A{r) 2 (apart from their 
parametrizations), the lens equation will depend on the 
metric coefficients S(r) and R(r) only. We will see below 
that many qualitative features of the lens equation are 
determined by the coefficient R(r) alone. 

For introducing our lens equation we have to fix two 
radius values ro and rs between r m ; n and r max . The 
index O stands for "observer", the index S stands for 
"source". We think of an observer at r — ro, <p = 0, 
■d = n/2. It is our goal to determine the appearance, on 
the observer's sky, of static light sources distributed on 
the sphere r = rs- 

By symmetry, we may restrict to the plane = n/2. 
We consider past-oriented lightlike geodesies that start 
at time t = at the observer and terminate, at some 
time t = —T < which depends on the geodesic, some- 
where on the sphere r — rs- To each of those light rays 
we assign the angle 0, measured at the observer between 
the ray's tangent and the direction of d r , and the angle 
swept out by the azimuth coordinate along the ray 
on its way from the observer to the source, see Figure 
n The desired lens equation is an equation of the form 
!F(Q, $) = which relates image positions on the ob- 
server's sky, given by 0, to source positions in the space- 
time, given by $ modulo 2ir. We restrict to values be- 
tween — 7r and 7r; then |0| can be viewed as a colatitude 
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coordinate on the observer's celestial sphere. By symme- 
try, $) = must be equivalent to !F(—Q, -$) = 0. 
For a given angle 0, neither existence nor uniqueness of 
an angle $ with F(Q, $) = is guaranteed. Existence 
fails if the respective light ray never meets the sphere 
r = rs; uniqueness fails if it meets this sphere several 
times. In the latter case the observer sees two or more 
images of light sources at rs at the same point on the 
sky, one behind the other. We will refer to images which 
are covered by other images as to "hidden images" . The 
lens equation can be solved for 4>, thereby giving a lens 
map <— > <fr, only if hidden images do not exist (or are 
willfully ignored). 

To work out the lens equation we have to calculate the 
lightlike geodesies in the plane # = ir/2 of the metric 0, 
which is an elementary exercise. As a conformal factor 
has no influence on the lightlike geodesies (apart from 
their parametrization), they are solutions of the Eulcr- 
Lagrangc equations of the Lagrangian C = \(S(r) 2 r 2 + 
R(r) 2 ip 2 - i 2 ), i.e. 
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S(r)S'{r)r 2 + R{r)R'(r)£ 2 
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t = , 



0. 
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(4) 



where an overdot denotes differentiation with respect to 
the curve parameter s. As an aside we mention that, by 
0, a circular light ray exists at radius r p if and only if 
R'(r p ) = 0. Comparing this condition with the equiva- 
lent but less convenient eq. (33) in Atkinson's article 6] 
shows that it is advantageous to write the metric in the 
form • The relevance of circular light rays in view of 
lensing was discussed by Hasse and Perlick [23], also see 
Claudel, Virbhadra and Ellis [2l| for related results. 

To get the past-oriented light ray that starts at time 
t = at the observer in the direction determined by the 
angle we have to impose the initial conditions 
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For each 0, the initial value problem 0, ©, 0,0, 0, 
has a unique maximal solution 

r = r(0,s), (p = (p(&,s), t = -s (8) 

where s ranges from up to some s max (0). Every image 
on the oberver's sky of a light source at rs corresponds 
to a pair (0, $) such that 

?\s = r(0,T) and $ = 0(0,T) (9) 



with some parameter value T £ ] 0, s max (0) [ . In other 
words, we get the desired lens equation !F(Q, $) = if 
we eliminate T from the two equations 0. 

We get an explicit expression for the lens equation, and 
for the travel time T, by writing the functions r(0, s) and 
</>(0, s) in terms of integrals. From the constant of motion 

S(r) 2 r 2 + R(r) 2 <p 2 - i 2 = (10) 

we find, with the help of 0, 0,0, 0, 

S{r) 2 R(r) 2 f 2 = R(r) 2 - i?(r o ) 2 sin 2 . (11) 

If f does not change sign, integration of 111|) yields 

|cos©| r {e - s) R(r)S{r)dr 



cos 



R{r) 2 - i?(r o ) 2 sin 2 



(12) 



With r(0, s) known, (p — 0(0, s) is determined by inte- 
grating with 0, 



0(0, a) 



R{ro) sin ds 
o i?(r(0,S)) 2 



(13) 



(|13f) can be rewritten as an integral over r, with r — 
dr/ds substituted from (|ll|l . This gives us the lens equa- 
tion in the form 

= |cos0| f rs R{r ) sin QS{r)dr 
cos© ./ 



R(r) \JR(r) 2 - R(r ) 2 sm 2 Q 

If f changes sign, (|T2"|> has to be replaced by a piece- 
wise integration. Similarly, the substitution from the 
s-integration in Ijl3(l to an r-integration must be done 
piecewise. In this case, the lens equation is not of the 
form l|14|) ; in particular, it is not guaranteed that the 
lens equation can be solved for $. In any case, we get 
exact integral expressions for the lens equation, and for 
the travel time T, from which all relevant lensing features 
can be determined in a way that is not less explicit than 
the quasi-Newtonian approximation formalism. This will 
be demonstrated by two examples in Section In Sub- 
section we treat a particularly simple example where 
the metric coefficients R(r) and S(r) are analytic and 
the integral Ijl2|l can be explicitly calculated in terms of 
elementary functions. In this case it suffices to calcu- 
late <|12f) for arbitrarily small s to get the whole function 
r(0, s) by maximal analytic extension; i.e., in this case it 
is not necessary to determine the points where r changes 
sign and to perform a piecewise integration. 



III. DISCUSSION OF THE LENS EQUATION 

In the first part of this section we want to discuss for 
which values of the lens equation <&) = admits 
a solution. In other words, we want to determine which 
part of the observer's sky is covered by the light sources 
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= TT ~ — 7T- 




e = o 



e = -tt/2 



FIG. 2: In the case ro < rs, the light sources cover the non- 
shaded region on the observer's sky, with 8i and Sn given 
by 11511 and 1161 1. respectively. (For rs < ro one gets an 
analogous picture, with O = and © = tt ~ — tt interchang- 
ing their roles.) If the spacetime is asymptotically flat and 
if rs is sufficiently large, this picture gives the socalled "es- 
cape cones" which have been calculated by Synge |22| for the 
Schwarzschild spacetime and by Pande and Durgapal 23] for 
an unspecified asymptotically flat spherically symmetric and 
static spacetime. 



distributed at r = r$- We restrict to the case ro < rs- 
(The results for the case ro > rs follow immediately 
from our discussion; we just have to make a coordinate 
transformation r — ► — r and, correspondingly, to change 
into 7T — 0. The case ro — rs can be treated by a limit 
procedure.) 

For a light ray with one end-point at ro and the other 
at rs the right-hand side of 111|) must be non-negative 
for all r between ro and rs- This condition restricts the 
possible values of O by sin 2 6 < sin 2 Si where 



sin Si = inf { R(r)/R(ro) \ ro < r < rs } 



(15) 



Note that our assumptions guarantee that this infimum 
is strictly positiv, < Si < 7t/2. 

Furthermore, a light ray with tt/2 < |0| < tt can arrive 
at rs only if it passes through a minimal radius value 
p(0) < ro- As JTTJl requires i?(p(0)) 2 = ,R(r o ) 2 sin 2 0, 
this can be true only if sin 2 > sin 2 Sn where 

sin Sn = inf { R(r)/R(r ) \ r min <r <r } , (16) 

< Si i < tt/2. So in general the light sources at rs cover 
on the observer's sky a disk of angular radius Si around 
the pole = and, if Sn < 5i, in addition a ring of 
angular width Si — Sn around the pole = tt ~ — tt, see 
FigureEl The two domains join if Si — tt/2. 

We see that the allowed values of are determined 
by the metric coefficient R alone. We will now demon- 
strate that R alone also determines the occurrence or 



non-occurrence of hidden images. Hidden images occur 
if a light ray from ro intersects the sphere r = rs at least 
two times; between these two intersections it must pass 
through a maximal radius tr(0) > rs which, by (|11|) . has 

to satisfy R(cr(Q)) = R(r ) 2 sin 0. Such a radius u(0) 
exists for all with sin 2 © > C 2 where 



C = inf { R{r)/R(r ) \r s <r < r max } . 



(17) 



As is restricted by sin 2 © < sin 2 <5/, hidden images can- 
not occur if C > sin Si. The latter condition is satisfied 
in asymptotically flat spacetimes, where R(r) — > oo for 
r — > r max (= oo), if we choose rs sufficiently large. This 
is the reason why in the more special situation of the 
Virbhadra-Ellis lens equation ^2] hidden images cannot 
occur. 

In the rest of this section we discuss the question of 
multiple imaging and the occurrence of Einstein rings. 
For a light source at r — rs, f — (fo, "d — tt/2 with 
< \<pq\ < tt, images on the observer's sky are in one-to- 
one correspondence with solutions © of the equation 



F{®,<Pq + 2mr) =0 



(18) 



with n € Z. We call the integer n the "winding number" 
of the corresponding light ray. An image with n = 
is called "primary" and an image with n — —fa/\fo\ is 
called "secondary" . Images with other values of n corre- 
spond to light rays that make at least one full turn and 
have been termed "relativistic" by Virbhadra and Ellis 
[T^ . Note that different images of a light source may 
have the same winding number. 

If we send ipo to or to tt ~ —tt, solutions of equation 
(|18|l with sin© ^ come in pairs (0,-0). By spherical 
symmetry, every such pair gives rise to an Einstein ring. 
There are as many Einstein rings as the equation 



He, 



= 



(19) 



admits solutions with positive integers i. Even integers i 
correspond to Einstein rings of the source at ipo = 0, and 
odd integers i correspond to Einstein rings of the source 

at ifo — TT ~ — TT. 



IV. OBSERVABLES 

To each solution (0, $) of the lens equation we can as- 
sign redshift, travel time, apparent brightness and image 
distortion. 



A. Redshift 

The general redshift formula for static metrics (see, 
e.g., Straumann [24]], P- 97) specified to metrics of the 
form ||TJ says that the redshift z is given by 



1 + z 



A(r s ) 



(20) 
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if the observer's worldline is a £-line at r = ro and the 
source's worldline is a i-line at r — rg. In our situation 
ro and rg are fixed, so the redshift is a constant. 



B. Travel time 

Recall that (9, $) is a solution of the lens equation if 
and only if there is a parameter T such that the equations 
@ hold. This assigns a travel time T to each solution 
(0, $) of the lens equation. If there are no hidden im- 
ages, the equation r(9, T) = rg gives T as a single- valued 
function of 0. 



C. Angular diameter distance 

Quite generally, determination of the angular diame- 
ter distance requires solving the Sachs equations for the 
optical scalars along lightlike geodesies, sec e.g. Schnei- 
der, Ehlers and Falco For the Schwarzschild met- 
ric, this has been explicitly worked out by Dwivedi and 
Kantowski [2f| . Their method easily carries over to arbi- 
trary spherically symmetric and static spacetimes as was 
demonstrated by Dyer [26j . In what follows we give a re- 
formulation of these results in terms of our lens equation. 

To that end we fix a solution (0, $) of the lens equation 
and thereby a (past-oriented) light ray from the observer 
at ro to a light source at rg. Around this ray, we consider 
an infinitcsimally thin bundle of neighboring rays, with 
vertex at the observer. The angular diameter distance 
is defined as the square-root of the ratio between the 
cross-sectional area of this bundle at the light source and 
the opening solid angle at the observer. Owing to the 
symmetry of our situation there are two preferred spatial 
directions perpendicular to the ray: a radial direction 
(along a meridian on the observer's sky) and a tangential 
direction (along a circle of equal latitude on the observer's 
sky). Therefore, the angular diameter distance naturally 
comes about as a product of a radial part and a tangential 
part. 

To calculate the radial part, we consider the infinitesi- 
mally neighboring ray which corresponds to an infinitcs- 
imally neighboring solution (0 + dQ, $ + d$) of the lens 
equation, i.e. g?<I> and dQ satisfy 



dQ 



(0,$) 



d$ { ' J dQ 



= 



We define the radial angular diameter distance as 

d$ 



A{rg)R(r s ) cos ip 



dQ 



(21) 



(22) 



with ip given by Figure^ i- e -> A(rg) R(rg) cos ip d& mea- 
sures, in the direction perpendicular to the original ray, 
how far the neighboring ray is away. By © , © and , 
ip must satisfy 



With ip given by ((23 and d$/dO given by J2TJ|, D£ ng is 
determined by \2'2\) for every solution (0,$) of the lens 
equation. Note that -C^ng is singular at those solutions 
of the lens equation where dT jd<& has a zero. If the lens 
equation can be solved for $, we may view -Dang as a 
(single- valued) function of 0. 

To calculate the tangential part we consider an in- 
finitesimally neighboring light ray that results by ap- 
plying a rotation around the axis tp = 0, d — tt/2. 
Such rotations are generated by the Killing vector field 
K = simpdv + cot •& cos ip d v . At points with — 
tt/2, this Killing vector field takes the form K — 
sinipA(r) R(r) g(d#, d$)~ 1 / 2 d§. Hence, if we rotate by 
an infinitesimal angle df3, the neighboring ray intersects 
the sphere r = rg at a distance A(rg) R(rg) sin <I> d{3 
from the original ray. Relating this distance to the angle 
sin d(3 between the two rays at the observer gives the 
tangential angular diameter distance 



D, 



= A(rg)R{r s )- 



sin<I> 



sin0 ■ ^ 
By this equation, -D ang is uniquely determined for each 
solution (0, $) of the lens equation. Again, -D ang may 
be viewed as a function of if the lens equation can be 
solved for <I>. 

D^ng an d ^Lig together give the (averaged) angular 
diameter distance or area distance 



i'a.ntr — \l I -Dang ^ang 



(25) 



Note that both -D^ng an d -D^ng ma y be negative. Images 
with -Da, lg -^ang > are said to have even parity and 
images with Da ng D* ng < are said to have odd parity. 
Images with odd parity show the neighborhood of the 
light source side-inverted in comparison to images with 
even parity. 

A solution (0, $) of the lens equation is called a ra- 
dial critical point if D 7 &ng = and a tangential critical 
point if -D^ ng = 0. The latter condition is equivalent to 
sin$ = and sin© ^ 0, i.e., to the occurrence of an 
Einstein ring. Note that (radial and tangential) critical 
points come in pairs, (0,$) and (—0,—$). Every such 
pair corresponds to a circle of equal latitude on the ob- 
server's sky which may be called a (radial or tangential) 
critical circle, as in the quasi-Newtonian approximation 
formalism, see Schneider, Ehlers and Falco Q], p. 233. 
In the quasi-Newtonian formalism one usually introduces 
the inverse magnification factors 1/ n r and as substi- 
tutes for DL,„ and -D* . In our situation, where there 
is no flat background metric, not even asymptotically, 
the magnification factors cannot be defined in a reason- 
able way, but working with -Da ng and D^ng is completely 
satisfactory. 



D. Luminosity distance 



-R(ro)sin0 = R(rg) sin ip . 



In arbitrary spacetimes, the angular diameter distance 
(23) £*ang is related to the (uncorrected) luminosity distance 
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Di um by the universal formula Z?i um = (1 + z) 2 -D all g) see, 
e.g., Schneider, Ehlers and Falco , eq. (3.80). With 
the redshift z given by l|20[l and the angular diameter 
distance D ang given by we find 



A ( r o) 2 hjyr D t 
A(r s ) 2 V 1 ans ang 



(26) 



For an isotropically radiating light source with bolo- 
metric luminosity L, the total flux at the observer is 
L/(47r_Dj 2 um ), see again Schneider, Ehlers and Falco 0, 
eq. (3.79). Hence, if we distribute standard candles at 
r = rg, their apparent brightness on the observer's sky 
is proportional to -Di~ ro - 



E. Image distortion 



D ln S and D* ng 



immediately give the apparent distor- 
tion of images. For the sake of illustration, we may 
think of small spheres, with infinitesimal diameter dD, 
distributed with their centers at r — r$ . By definition of 
D^ ng and , each solution (<d, $) of the lens equation 
corresponds to an elliptical image of such a sphere on 
the observer's sky, with the radial (meridional) diameter 
of the ellipse equal to |d8| = \dD/D^ ng \ and with the 
tangential (latitudinal) diameter of the ellipse equal to 
\d(3\ = \dD I '-Dangl- Thus, we may use the ellipticity 

\de\ - \dp\ \Dl ng 



\de\ 



\Dl 



(27) 



as a measure for image distortion. 



EXAMPLES 



7T J k 

.'J;; 










2n- 




















-7T -7r/2^--^" 


n/2 7T 


/ „jp 
















/ :S;i 




/ —n/k 







e 



FIG. 3: The lens equation Il.'-i2fl for the Barriola-Vilenkin 
monopole with k — 1/3.7 and ro = 0.77 rs. Intersections 
with the dashed lines indicate Einstein rings. Shading distin- 
guishes sources with N(k) + 1 — 4 images from sources with 
N(k) = 3 images. 
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A. Lensing by a Barriola-Vilenkin monopole 

We consider the metric 

g = dr 2 + k 2 r 2 {d$ 2 + sin 2 §dip 2 ) - dt 2 (28) 

where k is a positive constant. For k = 1, this is just 
Minkowski spacetime in spherical coordinates. For k < 1 , 
there is a deficit solid angle and a singularity at r = ; 
the plane t = const., i? = 7r/2 has the geometry of a 
cone. Similarly, for k > 1 there is a surplus solid angle 
and a singularity at r = 0. For k ^ 1, the metric is 
non-flat. The Einstein tensor has non- vanishing compo- 
nents Gtt = —G rr — (1 — fc 2 )/r 2 , so the weak energy 
condition is satisfied (without a cosmological constant) 
if and only if k < 1. In that case it was shown by Bar- 
riola and Vilenkin 27] that the metric may be viewed 
as a model for the spacetime around a monopole result- 
ing from breaking a global 0(3) symmetry. To within 
the weak-field approximation, basic features of lensing 
by such a monopole were discussed in the original paper 
by Barriola and Vilenkin [27j and also by Durrer (2g • In 
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/ 3;; 




' —n/k 







FIG. 4: The lens equation (1321 for the Barriola-Vilenkin 
monopole with k — 1/3.7 and rs = 0.77 ro- Other than 
in the case ro < rs, the lens equation cannot be solved for 
<&, i.e., there are hidden images. The angle a is determined 
by rs = rosin a and n/2 < a < n. The picture shows 
that all Einstein rings, indicated by intersections with the 
dashed lines, are hidden. Shading distinguishes sources with 
N(k) + 1 — 4 images from sources with N(k) — 3 images; only 
one image of each source is non-hidden. 



what follows we give a detailed account in terms of our 
exact len s eq uation. Note that the Virbhadra-EUis lens 
equation [T^] is not applicable to this case, at least not 
without modification, because for k ^ 1 the spacetime is 
not asymptotically flat in the usual sense. 

Comparison of l|28|) with shows that the metric co- 
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D r D 



rs + r a 




e 



FIG. 5: Tangential angular diameter distance Dl ng (dashed) 
and radial angular diameter distance -D ang (solid) as functions 
of Q for the Barriola-Vilenkin monopole with k — 1/3.7 and 
tq — 0.77 rs- The dashed curve also gives the travel time 
T — -Ding- Each zero of Dang indicates an Einstein ring. 
Where D ang and -D^ng have the same sign, the images have 
even parity; where -D ang and Dl ug have different signs, the 
images have odd parity, i.e., the neighborhood of each light 
source is shown side-inverted. 



efficients are given by 

A(r) = 1 , S{r) = 1 , R(r) = kr (29) 



on the interval 



= < r < r n 



= oo. With these 



metric coefficients, the integrals 11121 and (1131 can be 
calculated in an elementary fashion, yielding the solution 
to the initial value problem in the form 



r(6,s) 



+ 2 tq s cos + s 2 



(30) 



ro 



! COS 



1 sin9 / 

<p(0,s) = — — ■ arccos — , 

k |sm0| V r o + 2r scose + s 2 

(31) 

For sin0 ^ 0, s ranges from to oo, so \ip\ — \<fi(0,s)\ 
ranges from to |0|/fc. Eliminating T from the two 
equations © gives the lens equation, 

r s sin(0 - fc$) - r sin© = , (32) 

which is to be considered on the domain 



- 7T < < 7T , 



fc|$| < |0| 



(33) 



We restrict to the case that the integer N(k) defined by 
N(k) < 1/k < N(k) + 1 is odd. The lens equation is 
plotted for the case ro < rs in Figure OH and for the 
case ro > r§ in Figure (For producing the pictures 
we have chosen k such that N(k) — 3.) In either case 
we find that there are N(k) Einstein rings. For a light 




+ m 



e 



FIG. 6: Apparent brightness of standard candles for the 
Barriola-Vilenkin monopole with k — 1/3.7 and ro = 0.77 rs- 
Instead of D^^, which is proportional to the energy flux 
from point sources, we have plotted the magnitude m = 
2.5 log 10 (_D! 2 um ) + mo used by astronomners. The constant 
mo has been chosen such that m = at O = 0. Einstein rings 
are infinitely bright in the ray optical treatment, m = — oo. 



source at r — rs,ip = <^0:^ = f/2 with < \fo\ < tt 
there are N(k) images if \(po\ < (l + N(k) — 1/fcW (non- 
shaded regions in Figures |21 and and N(k) + 1 images 
otherwise (shaded regions in Figures and 0J . 

In the following we concentrate on the case ro < rs- 
Then the lens equation can be solved for giving a lens 
map 







-J 



— arcsinf — sin© 
V r s 



(34) 



on the domain —tt < < 7r, i.e., in the notation of 
Figure |21 we have Si = n/2 and 8u = 0. In the case 
k = 1 (flat spacetime) , the lens map can be continuously 
extended into the point — 7r ~ — tt (and by periodicity 
onto the full circle, € M mod 2tt). For k ^ 1 the ray 
with = n ~ — tt cannot pass through the singularity 
at r — 0, so the lens map is not to be extended beyond 
the open interval ] — tt, tt [ . On this interval, $ increases 
monotonously from —n/k to n/k, see Figure|3] Thus, the 
range of the lens map is independent of ro and rs- The 
ratio ro /rs influences the shape of the graph of the lens 
map in the following way. For ro /rs — * it becomes 
a straight line, $ = Q/k. For ro/rs — ► 1 it becomes 
a broken straight line, $ = for |0| < 7r/2 and $ = 
(20 - ir)/k for |0 - tt\ < tt/2. Note that linearity of 
the lens map implies that the angular distance on the 
observer's sky between consecutive images is the same 
for all light sources. 

As $ is a (single- valued) function of 0, so are all ob- 
servables. By evaluating the formulas derived in Section 
IV] for the case at hand we find 



2 = 0, 



(35) 



T 



Ding = r O COS + 



^sin 2 © , (36) 
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FIG. 7: The ellipticity e for the Barriola-Vilenkin monopolc 
with k — 1/3.7 and ro = 0.77 rs- If small spheres are dis- 
tributed with their centers at r — rs, this function describes 
their distortion into ellipses on the observer's sky. e — 
indicates circular images, e = — 1 images with no tangential 
extension (radial critical points) , and e = 1 images with no ra- 
dial extension (tangential critical points, i.e. Einstein rings). 



D* 

ang 



kr s 



sin I \ ( — arcsin( 
sin 6 \ k \ 



ro sin 

rs 



(37) 

which gives D\ um and £ as functions of via 1(26(1 and 
1)27(1. The observables are plotted in Figures [3 El an d 13 



B. Lensing by an Ellis wormhole 

The metric 




(-) 



FIG. 8: The lens equation g2J for the Ellis wormhole. The 
angle <5j(< 7r/2) is defined in 1401 . There are infinitely many 
Einstein rings, indicated by intersections with the dashed 
lines. 



e 



g = dr 2 + (r 2 + a 2 ) {d$ 2 + sin 2 $ dtp 2 ) - dt 2 , (38) 

where a is a positive constant, is an example for a 
traversable wormhole of the Morris-Thorne class, see Box 
2 in Morris and Thorne l29l . It was investigated, already 
in the 1970s, by Ellis [3(| who called it a "drainhole". 
Lensing in the Ellis spacetime was discussed, in a scat- 
tering formalism assuming that observer and light source 
are at infinity, by Chetouani and Clement |3l|. In the 
following we give a detailed account of lensing in this 
spacetime with the help of our lens equation. 
By comparison of l|38|l with we find 

A{r) = 1 , S(r) = 1 , R(r) = \Jr 2 + a 2 . 

(39) 

where the radius coordinate r ranges from r m j n = — oo 
to r max = oo. (We do not identify the region where r is 
positive with the region where r is negative.) The func- 
tion R has a minimum at r — 0, thereby indicating the 
existence of circular light rays at the neck of the worm- 
hole. In the following we consider the case that observer 
and light sources are on different sides of the neck of the 
wormhole, — oo < ro < < rs < oo. 



FIG. 9: The travel time T as a function of Q for the Ellis 
wormhole. The axis O = is met at To = rs — ro, see 1441 . 

As a first step, we determine for which angles the 
lens equation admits a solution, recall Figure [21 In the 
case at hand, the angles 8i and Su defined by (|15|) and 
are given by sin 8j = a/ \J r 2 Q + a 2 and Su — n/2 > 
Si. Hence, the lens equation admits a solution for all 
angles with |0| < 7r/2 and 

a 2 

sin 2 © < sin 2 5/ = -„ . (40) 

r 2 Q + a 2 

Light sources distributed at r = rs illuminate a disk 
of angular radius S] < ir/2 on the observer's sky. The 
apparent rim of the disk corresponds to light rays that 
spiral asymptotically towards r = 0. As the constant C 
defined by |(T7Jl satisfies C = y/ (r| + a 2 ) /(r + a 2 ) > 
a/y/rfy + a 2 = sin<5j, there are no hidden images, i.e., 
the lens equation can be solved for <&. 

With restricted by (|4"0"|l we read from (fTT|l that f 
has no zeros along a ray that starts at r — ro and passes 
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J -'ang ■> ^ang 



D 




FIG. 10: -D^ng (dashed) and -D ang (solid) as functions of O 
for the Ellis wormhole. At each Einstein ring -D an g has a 
zero. To make the oscillatory behaviour of -D* visible, we 



use x — — log(l — |0|/5rJ instead of O on the abscissa; |0| 
ranges from to 5i if x ranges from to oo. The axis is met at 
Do = \/ r s + a2 \/ r o + 0,2 a ~ 1 ^arctan(rs/a) — arctan(ro/a)^ , 
see (|i5l and ijlfTjl . 

m = 2.5log 10 (Df um ) + m 



30 



20 
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FIG. 11: The apparent brightness, measured in terms of mag- 
nitude m = 2.5 log 10 (-D[ 2 um ) + mo, as a function of for the 
Ellis wormhole. The constant mo has been chosen such that 
m = at = 0. We use the same coordinate x as in Figure 
HOI on the abscissa. At each Einstein ring there is a singularity, 



through r = r$- Hence, ()12t gives us r(0, s) and thereby 
the travel time in terms of an elliptic integral, 



■(e, s ) 



\Jr 2 + a 2 dr 



a 2 cos 2 9 — rQsin 2 



(41) 



Similarly, (|14l) gives us $ as a (single-valued) function 
of in terms of an elliptic integral and thereby the lens 
equation, 



\/tq + a 2 sin O dr 



Vr 1 



a 2 cos 2 



r 2 Q sin 2 6 



see Figure |S1 As $ increases monotonously from — oo to 
oo on the domain ] — 8i, Si [, there are infinitely many 
Einstein rings whose angular radii converge to Si. If we 
fix a light source at <p = ipo with < \<po\ < tt, Figure 
|S] gives us infinitely many images which can be charac- 
terized in the following way. For every n £ Z there is a 
unique 0„ G ] — Sj, Sj [ such that T(Q n , tpa + 2nn) = 0, 
and 0„ — > ± Si for n — » ± oo. 

As $ is a (single- valued) function of 0, so are all ob- 
servables. By evaluating the formulas derived in Section 
IIVI for the case at hand we find 



0. 



(43) 



\/ r 2 + a 2 dr 



a 2 cos 2 



r 2 sin 2 



(44) 



Dl 



a 2 cos 2 — r^sin 2 © y/ r 2 Q + a 2 V r 2 + a 2 cos dr 



(45) 



a 2 cos 2 



r 2 Q sin 2 © 
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D* 

ang 



sin 9 



J rr 



\J r Q + a 2 sin 9 dr 



V r 2 + , 



a 2 cos 2 9 



r 2 Q sin 9 



(46) 



which gives D\ um and e as functions of 9 via and 
(|27j) . The observables are plotted in Figures 131 ITU1 ITT1 
and IHfl 

As there are infinitely many Einstein rings whose angu- 
lar radii converge to Si, the tangential angular diameter 
distance must have infinitely many zeros that converge to 
Si. This is difficult to show in a picture unless one trans- 
forms 9 into a new coordinate x that goes to infinity 
for 9 — > Si. If one chooses a logarithmic transformation 
formula, as has been done in Figure one sees that 
in terms of the new coordinate x the Einstein rings be- 
come equidistant. This feature is not particular to the 
Ellis wormhole; <!> as a function of 9 always diverges log- 
arithmically when a circular light ray at a radius r p with 
R'{r p ) = and R"{r p ]> is approached. The proof can 
be taken from Bozza |l8j . 

One may also treat the case that observer and light 
sources are on the same side of the neck of the wormhole. 
If the observer is closer to the neck than the light sources, 
< ro < rs < oo or — oo < rs < ro < 0, the results 
are quite similar to the case above. The only difference 
is in the fact that the light sources appear as a disk of 
radius bigger than 7r/2, i.e., the disk covers more than 
one hemisphere. If the observer is farther from the neck 
than the light sources, < rs < ro < oo or — oo < ro < 
rs < 0, there are hidden images, i.e., one does not get a 
single- valued lens map 6 n $. 

The qualitative features of lcnsing by an Ellis worm- 
hole are very similar to the qualitative features of lensing 
by a Schwarzschild black hole. The radii r m i n = — oo, 
r p = 0, r max = oo in the Ellis case correspond respec- 
tively to the radii r m j n = 2m, r p = 3m, r max = oo 
in the Schwarzschild case. As a matter of fact, we en- 
counter these same features whenever the function R has 
one minimum, R'{r p ) — and R"{r p ) > 0, and no other 
cxtrema on the considered interval. 



VI. CONCLUDING REMARKS 

The lens equation and the formulas for redshift, travel 
time and radial angular diameter distance used in this pa- 
per refer to lightlike geodesies of the (2 + l)-dimensional 
metric A(r) 2 ( S (r) 2 dr 2 + R(r) 2 dip 2 — dt 2 ), indepen- 



dently of whether this metric results from restricting a 
spherically symmetric and static spacetime to the equa- 
torial plane. Therefore, these results apply equally well 
to the plane z = const, of a cylindrically symmetric and 
static spacetime and, of course, to genuinely (2 + 1)- 
dimensional spacetimes with the assumed symmetries 
such as the BTZ black hole. (For lightlike - and time- 
like - geodesies in the metric of the BTZ black hole 



FIG. 12: The ellipticity e as a function of for the Ellis 
wormhole. Again we have chosen the coordinate x from Fig- 
ure ED on the abscissa, e is equal to 1 at the Einstein rings 
and smaller than 1 between the Einstein rings. However, be- 
yond the third Einstein ring e stays so close to 1 that in the 
picture it looks like a constant function. 



see Cruz, Martinez and Pena 32].) E.g., the metric 
dr 2 + k 2 r 2 dip 2 — dt 2 results not only by restricting the 
spacetime of a Barriola-Vilenkin monopole to the plane 
■d = 7r/2, as discussed in Subsection IV Al but also by 
restricting the cylindrically symmetric and static metric 
dr 2 + k 2 r 2 dip 2 +dz 2 - dt 2 to the plane = 0. The latter 
metric is well-known to describe the spacetime around a 
static string, see Vilenkin |33] , Gott [34J and Hiscock 
and was investigated in detail already by Marder [3£ 
Hence, if re-interpreted appropriately, the results of Sub- 
section lY Al applv to light rays in the plane perpendicular 
to a static string. For treating all light rays in a cylindri- 
cally symmetric and static spacetime one may introduce 
a modified lens equation, replacing the sphere r = rs 
with a cylinder. 
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